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Operation of a Z Pinch 
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Heating Mechanisms Included in Model 
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Fission Power Dominates as 
Neutron Count 
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• Parameter space for 
ignition 

• Greatly broadened with 
embedded magnetic field 

• Marginally improved with 
6 Li and thorium liners 

• Significantly enhanced 
with uranium liners ( 235 U 
and 23 8 U) 






Parameters within net Power Increase 
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♦ Choose p R = 10 -5 at 15 keV. 

♦ Let R = 1 mm, thickness of uranium liner is 5 mm, length of target is 2 cm 

♦ Density of DT target is 0.1 kg/m 3 

♦ Total energy in DT target is only ~5 kJ, 1% of Charger 1 stored energy 
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Our Approach: Solve Maxwell's Equations Coupled to 
Multifluid (Ions, Electrons, Neutrals) Equations of Motion 



Maxwell's Equations 


• Solve with Smooth Particle Electromagnetic Variant of Finite-Difference 
Time Domain (FDTD) method 

• FDTD well documented, highly accurate grid-based method for analyzing the 
time evolution of electric and magnetic fields, utilized in PIC codes 

• Can interpolate charged fluid particles to grid to model conductivity or charge 
and current density 


Multifluid Equations of Motions 


• Solve with Smooth Particle Hydrodynamics (SPH) 

• Gridless Lagrangian technique 

• Vacuum/plasma boundary well defined 

• Leverage same engine as Maxwell Equation Solver 
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Both methods yield to Vectorized' coding, making 
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multiprocessor (parallel) computing easy 








♦ Numerical method for approximating probability densities over a 
domain of particles. 

♦ Developed initially to model stars by Monaghan, et al, in 1977. 

♦ Currently used mostly in hydrodynamic modeling and CG effects in 
film and video games. 
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Integral interpolant: 

A\ (r) = J A(r')W(r — r',h)dr' 

A - quantity measured (density, 
temperature, etc.) 

W— differentiable kernel function 
dr' - volume differential 
h - smoothing length. 

Summation over mass elements: 



A j>j 

— W(r — r b , K) 
Pb 


Similar to density probability 
calculations. 

How quantities are accurately calculated 
with a small particle domain. 
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Equations of motion (completed) 
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Transport effects, 
which can be based 
on nonequilibrium 
distribution functions 
(kappa and power law) 
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Initial Pulsed Nozzle Model 
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Test thermal expansion of 
gas nozzle with various 
initial conditions 

• Nozzle geometry 

• Gas 


• Temperature 

• Density 

• Radius 

(m) 

• Length 

• Composition 

Lays ground work and 
expectations for magnetic 
nozzle 


Preliminary results 
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NIAC Phase I Goals 


Crewed Mars Mission Concept 



D eutei i u in -T i iti u in Tan k 
(2.4 in did.) (4 pics) 



Stacked Capacitor Module (2) 
(10 in long x 3.0 in x 7.2 mi 
(8 pics) 
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Radiation Shield 
(.25 in thick) 
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(1.8 m x 1.8 m x 2.4 m) 
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Trans hah 


Surface 


Two-Sided Med. Temp. Radiators 
(008 m2 total area) 


Lander 


Habitat 


Two-Sided High Temp. Radiators 
(1910 m2 total area) 


4-Pod RCS 700 Ibf MR 80B 3 Hydrazine 
Thrusters, RCS Tank (.80 in dia.) and RCS 
Helium Pressurant Tank (.04 in dia.) 

(8 pics - 4 Aft and 4 Forward) 
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Mission Concepts 
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Figure 3 Mars 90 Day Transfer Trajectories 


» Engine 

• Isp = 19,400 sec 

• T = 38kN 

• 10 Hz pulse freq. 

* Vehicle 

• M dry = 552 mT 

• M pay = 150 mT 

• 30% MGA 

Polsgrove, T. et al. Design of Z-Pinch and Dense 
Plasma Focus Powered Vehicles, 2010 AIAA 
Aerospace Sciences Meeting 


22 




















Mating SPFMaX and MCNP 


♦ SPFMax gives 

• Ability to model 3d effects 

• Can propagate magnetic fields in vacuum 

• Easily editable 

♦ MCNP 

• Track neutron life, fission reactions 

• Flexible geometries 

♦ Second half of NIAC is to run codes concurrently 

• synchronize neutron population vs. time 

• Optimize energy output 

- As function of geometry 

- As function of composition 

- Mix of UF6, D-T 

- Lithium liner thicknesses 
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Single turn Magnetic Nozzle 
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Gasdynamic nozzle 
performance to be compared 
with magnetic nozzle to assess 
loss mechanisms in magnetic 
nozzles, e.g. 

• Field/plasma instabilities 

• Plasma detachment 




NIAC Phase II Experimental Options 


Charger - 1 


♦ A test facility for high power and 
thermonuclear fusion propulsion 
concepts, astrophysics modeling, 
radiation physics 

♦ Located in the UAH Aerophysics 
Lab at Redstone 

♦ The highest instantaneous pulsed 
power facility in academia - 572 
kJ (1 TW at 100 ns) 
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Experimental Roadpath 



♦ Methodology 

• Incremental improvements in experimental capability 

• Benchmark model with experimental data 

• Can also run any experiments below with lower power systems 

• Looking for comments and suggestions here! 
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Long Range Plans 



♦ Charger II 

• Construct breadboard PuFF system capable of 10-20 Hz operation 

- Upgrade to flight weight hardware - NASA 

- Optimize pulse for maximum power output - DOE 

- Astrodynamics, radiation protection, other research goals - Various 
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